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摘 要 采用 电场 对 微 液 滴 进 行 精确 操控 在 化 工 、 生 物 医学 、 能 源 和 环境 领域 具有 广阔 的 应 用 前 景 。 准 确 预测 外 加 电场 作用 
下 粘性 液 滴 的 变形 与 破裂 行为 ， 是 精确 操控 微 液 滴 行为 的 重要 前 提 。 本 文 基于 漏电 介质 模型 发 展 格子 Boltzmann 有 限 差 分 混 
合 数值 方法 对 恒 稳 电场 作用 下 的 液 滴 变形 与 破裂 行为 进行 研究 ,结果 表明 : 本 文 数值 方法 可 对 恒 稳 电场 作用 下 液 滴 的 小 变形 
进行 准确 预测 ， 验 证 了 该 模型 的 有 效 性 和 准确 性 ， 对 于 沿 垂直 于 电场 方向 变形 的 液 滴 ( 扁 平 型 ) ， 液 滴 随 电场 毛细 数 的 增加 呈 
现 四 种 不 同 的 形态 ， 椭 圆 形 、 类 椭圆 形 、 哑 铃 型 和 破裂 ; 对 于 沿 电场 方向 变形 的 液 滴 ( 扁 长 型 )， 其 变形 参数 随 电场 毛细 数 的 
变化 与 扁平 型 液 滴 显 著 不 同 ， 液 滴 呈 现 的 形态 包括 椭圆 形 、 类 椭圆 形 、 周 期 振荡 和 破裂 。 在 液 滴 变 形 阶段 ， 扁 平 或 忆 长 型 液 
滴 均 随 电 场 毛细 数 的 增加 而 变形 增 大 ; 液 滴 破 裂 阶段 ,破裂 所 需 时 间 均 随 电 场 毛细 数 的 增加 而 缩短 。 研 究 工 作为 深刻 认识 电 
场 作用 下 微 液 滴 变形 和 破裂 机 制 提供 理论 基础 。 
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Numerical Study of Microdroplet Deformation and Breakup under a Steady Electric Field 


CUI Yu-Tong, WANG Ning-Ning, LIU Hai-Hu 
(School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China) 


Abstract Manipulation of microdroplets in a precise manner by an external electric field has received much 
attention due to its importance in the fields of chemical engineering, biomedicine, energy and environment. In 
order to accurately manipulate the droplet behavior in microfluidic devices, it is necessary to study and gain an 
in-depth understanding of the deformation and breakup of a viscous droplet under an external electric field. In this 
paper, a hybrid method, in which the color lattice Boltzmann and the finite difference methods are coupled by the 
leaky dielectric model, is used to study the deformation and breakup of a droplet under a steady electric field. The 
hybrid method is first validated by simulating the droplet subject to a small deformation. We then investigate the 
influence of the intensity of the electric field, the ratio of the dielectric permittivity inside and outside of the 
droplet, as well as the ratio of the electric conductivity on the droplet behavior. It is found that the oblate droplet 
can exhibit four different shapes or states, i.e. ellipse, quasi-ellipse, dumbbell, and droplet breakup, which are 
strongly dependent on the electrical capillary number. Distinct from the oblate droplet, the prolate droplet does not 
exhibit the dumbbell shape but the periodic oscillation. For either prolate or oblate droplet, increasing electrical 
capillary number leads to an increased droplet deformation when the deformation occurs, and to a shorter time for 
the droplet to break up when the breakup occurs. This study provides the theoretical foundation for a deep 
understanding of droplet deformation and breakup in the presence of electric field. 
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0 引言 


微 液 滴 具 有 体积 小 、 低 扩散 、 无 交叉 污染 、 快 
速 的 反应 动力 学 、 高 通 量 的 分 析 潜 力 等 优点 。 随 着 
微 流 控 技术 的 不 断 进步 与 发 展 , 微 液 滴 技 术 在 化 学 、 
生物 医学 、 能 源 和 环境 等 领域 得 到 了 广泛 的 关注 和 
应 用 趾 。 例 如 ， 在 化 学 领域 ， 液 滴 微 流 控 正 逐步 取 
代 传 统 分 析 化 学 中 的 实验 室 ， 微 液 滴 越 来 越 多 地 作 
为 纳 升 化 学 的 反应 器 或 合成 器 ， 精 确 控制 反应 试剂 
的 分 子 浓度 和 反应 时 间 ， 从 而 自动 快速 地 完成 各 种 
复杂 的 化 学 反应 和 分 析 功 能 。 在 微 液 滴 的 应 用 中 ， 
液 滴 动 力学 行为 的 精确 操控 是 极其 重要 的 。 存 在 多 
种 方式 生成 和 操控 微 液 滴 ， 其 中 ， 电 场 力 操控 是 一 
种 有 效 的 非 接触 式 方式 ， 并 已 成 功 应 用 于 筋 化 、 喷 
墨 打印 、 增 强 液 滴 融 合 、 破 乳 等 行业 。 在 外 加 电场 
的 作用 下 , 悬浮 在 粘性 液体 中 的 液 滴 可 表现 出 变形 、 
运动 、 旋 转 和 破裂 站 等 复杂 行为 。 这 些 复杂 行为 主 
要 取决 于 电场 强度 和 流体 性 质 ， 包 括 粘性 、 界 面 张 
力 、 导 电 性 和 介 电 常数 等 。 因 此 ， 开 展 电场 作用 下 
液 滴 变形 与 破裂 行为 的 研究 既是 微 流 控 技 术 发 展 的 
内 在 要 求 ， 也 可 为 工业 应 用 电场 提供 理论 指导 和 技 
ASCH 

对 于 电场 作用 下 不 混 溶液 滴 变 形 运 动 的 研究 ， 
早期 开始 于 Taylor”! VAA Melcher 和 Taylor Zi 
和 实验 研究 工作 。Taylor 关于 液 滴 变 形 的 理论 研究 
基于 一 些 基本 假定 ， 即 液 滴 呈 中 性 、 准 静态 电场 、 
液 滴 与 运载 流体 均 为 漏电 介质 流体 、 液 滴 只 做 微小 
变形 以 及 液 滴 表 面 无 电荷 移动 ， 相 关 理 论 称 作 漏 电 
介质 模型 。 在 这 种 假设 条 件 下 , 液 滴 本 身 不 带电 荷 ， 
但 是 允许 电荷 在 液 滴 表 面 上 堆积 。 由 于 电荷 堆积 速 
度 比 流 体 的 运动 速度 快 得 多 ， 导 致 系统 电荷 出 现 伪 
定常 分 布 ， 电 场 满足 准 静 态 假设 。Hua “ACE 
漏电 介质 模型 和 理想 介质 模型 ， 采 用 前 向 追踪 法 对 
带电 液 滴 在 均匀 电场 中 的 变形 进行 了 研究 ， 证 实 该 
理论 可 对 小 变形 、 有 限 导 电 性 和 介 电 常数 之 差 的 液 
滴 在 恒 稳 电场 中 的 变形 做 出 有 效 的 预测 。 然 而 ， 
Torza 等 人 四 实验 研究 发 现 液 滴 变 形 较 大 时 ，Taylor 
模型 的 理论 解 与 实验 结果 存在 较 大 偏差 。 随 后 ， 
Ajayi! Xt Taylor 所 采用 的 线性 理论 进行 了 进一步 延 
伸 ， 提 出 了 含有 二 阶 项 的 理论 模型 ， 但 是 该 模型 仍 
未 从 根本 上 解决 理论 解 与 实验 值 之 间 偏 差 的 问题 。 

在 电场 的 影响 下 ， 液 滴 不 仅 产 生变 形 甚 至 可 能 
破裂 ;感应 电荷 在 液 滴 表面 聚集 ， 同 时 液 滴 内 部 和 
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外 部 的 粘性 流动 相互 耦合 ， 且 内 外 流体 之 间 存 在 界 
面 张 力 ， 对 理论 分 析 电 流体 动力 学 问题 提出 了 巨大 
的 挑战 。 随 着 计算 机 运行 速度 的 飞速 提高 ， 数 值 建 
模 和 模拟 已 成 为 研究 电场 作用 下 液 滴 动 力学 行为 的 
另 一 种 重要 手段 。Sherwood!" 通 过 边界 积分 法 求解 
控制 电场 的 拉 普 拉 斯 方程 和 控制 流 场 的 斯 托 克 斯 方 
程 研究 了 液 滴 在 电磁 场 内 的 破裂 。 基 于 漏电 介质 模 
型 ,Baygents 四 采用 类 似 的 方法 研究 了 均匀 电场 中 两 
液 滴 的 相互 作用 。 同 样 基 于 漏电 介质 模型 ，Feng 和 
Scott 采用 贾 辽 金 有 限 元 方法 建立 了 在 Stokes 流动 
条 件 和 有 限 雷 诺 数 下 液 滴 变 形 与 电场 强度 之 间 的 关 
系 。Lac All Homsy! 采用 边界 积分 法 对 三 维 液 滴 在 
恒 稳 电场 作用 下 的 变形 进行 了 数值 研究 。 Fernandez 
等 人 M1 分析 了 槽 道中 的 液 滴 在 电场 力作 用 下 的 变 
形 。Tomar 等 人 HM 使 用 Volume-of-Fluid 方法 模拟 了 
导电 液 滴 在 电场 作用 下 的 变形 行为 。 

近年 来 ， 格 子 Boltzmann 方法 (LBM) 已 发 展 成 
为 模拟 复杂 流体 流动 的 一 种 有 效 数 值 工具 ,LBM A 
有 天 生 的 并 行 特 性 ， 以 及 边界 处 理 简 单 、 程 序 易 于 
实施 等 优点 。 此 外 ,LBM 的 动 理论 特性 使 得 它 能 够 
方便 地 描述 不 同 相 间 的 相互 作用 , 自动 追踪 相 界 面 。 
因此 , LBM 特别 适合 于 模拟 微 液 滴 动 力学 行为 。 本 
文 基于 LBM 颜色 模型 ”由 对 两 相 流动 进行 求解 ， 
并 采用 有 限 差分 法 对 电场 控制 方程 进行 求解 ， 两 种 
方法 通过 漏电 介质 模型 进行 耦合 ， 对 恒 稳 电场 作用 
下 不 可 压 粘 性 液 滴 的 变形 与 破裂 行为 进行 数值 研 
究 ， 揭 示 电 场 强度 、 液 滴 内 外 介质 介 电 常数 比 和 电 
导 率 比 对 液 滴 行 为 的 影响 规律 。 
1 研究 对 象 及 方法 
1.1 研究 对 象 

将 初始 半径 为 尺 、 电 中 性 的 二 维 圆 形 液 滴 置 于 
两 平行 板 电 容器 所 形成 的 恒 稳 电场 下 ， 液 滴 位 于 两 
板 中 心 并 悬浮 于 运载 流体 中 。 根 据 对 称 性 ， 本 文 取 
液 滴水 平 轴线 以 上 为 计算 域 , 几何 模型 如 图 1 所 示 。 
上 极 板 具有 恒定 电势 @ ,计算 域 下 边界 为 水 平 对 称 
轴 ， 其 电势 为 零 ， 在 系统 中 存在 着 垂直 向 下 、 强 度 
为 E 的 恒 稳 电场 。 

液 滴 与 运载 流体 分 别 用 下 角 标 “让 及 “o” 表 示 ， 
其 物理 性 质 分 别 为 : 密度 P P: 电导 率 0;、0,; 
介 电 常数 n E: 以 及 动力 粘度 1、X, 。 通 过 它 
们 可 以 定义 以 下 无 量 纲 参数 ， 即 密度 比 
M=P/p,， 电 导 率 比 R=Go/0, ， 介 电 常 数 比 


QO=&/&,， 及 粘性 比 和 = /1 。 为 简单 起 见 ， 
本 文 假设 液 滴 与 运载 流体 具有 相同 的 密度 和 粘度 ， 
旦 忽略 重力 的 影响 。 液 滴 和 外 部 运载 流体 之 间 界 面 
张力 系数 为 y 。 


D 人 
(Po, Oo, Eo, klo) 
E YI 
e jely 
(Pi, i, €i, Li) R, 
NY 


PAN T 


a SR; _ 
图 1 恒 稳 电场 作用 下 液 滴 几何 模型 示意 图 (中 间 点 划 线 为 竖 
直方 向 对 称 轴 ) 


Fig.1 Sketch of the geometrical model of a droplet immersed in 
a second fluid and subject to a steady electric field(the middle 
dash-dot line represents the symmetry axis in the vertical 


direction) 


由 于 仅 在 液 滴 表 面 存 在 电荷 ， 所 以 电场 力 仅 作 
用 在 液 滴 表 面 ， 并 通过 电场 毛细 数 Car KWE RHE 
度 与 界面 张力 的 相对 大 小 ， 其 定义 为 : 


_EsR 
4 


Taylor 采用 线性 模型 中 计算 漏电 介质 液 滴 在 恒 
稳 电 场 中 的 稳定 形态 ， 发 现 液 滴 在 电场 作用 下 会 出 
现 扁平 型 或 局 长 型 。 根据 Taylor 理论 ， 引 入 变形 参 
ND 定量 描述 液 滴 的 变形 : 

2 L-B 2 
L+B 

其 中 , LAB 分别 为 液 滴 的 长 轴 和 短 轴 长 度 。 电 场 
作用 下 液 滴 的 变形 主要 取决 于 以 下 无 量 纲 参数 : 电 
导 率 之 比 R， 介 电 常 数 之 比 O 和 电场 毛细 数 Cap. 

针对 恒 稳 电场 作用 下 二 维 液 滴 的 变形 运动 ， 
Feng 等 人 借鉴 Taylor 模 型 推导 出 变形 参数 的 理论 解 
D5 为 


Ca, (1) 


P 
7/ 30+RyP F 


其 中 f =R? +R+1-30 为 指示 函数 。 可 根据 指示 
函数 的 正 负 判 定 液 滴 的 变形 方式 ，_f > OM, OR 


(3) 
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沿 电场 强度 方向 变形 ， 为 扁 长 型 ，f4 <0 时 ， 液 滴 
治 垂直 于 电场 强度 的 方向 变形 ， 为 书 平 型 。 

值得 指出 ， 计 算 区 域 的 大 小 取决 于 液 滴 的 变形 
情况 ， 设 置 原则 为 在 尽量 节省 计算 成 本 的 前 提 下 ， 
液 滴 两 端 与 计算 区 域 边界 始终 保持 合理 的 距离 。 在 
无 特殊 说 明 的 情况 下 ， 本 文 对 称 计算 区 域 的 尺寸 给 
定 为 8R, x4R,。 
1.2 研究 方法 


本 文采 用 格子 Boltzmann 颜色 模型 "对 两 相 流 
动 进 行 求解 。 在 颜色 模型 中 ， 采 用 不 同 的 颜色 来 区 
分 不 同 的 流体 ， 不 同 流体 之 间 的 相互 作用 通过 引入 
颜色 梯度 来 实现 , 采用 Brackbill 等 人 59 的 连续 表面 
力 模 型 建 模 界面 张力 , 并 通过 Guo 等 上 的 体积 力 格 
式 引 入 到 模型 中 。 然 而 ， 该 处 理 并 未 保证 两 相 流 体 
互 不 相 溶 的 特性 。 为 促使 液 滴 与 运载 流体 两 相 分 离 
并 维持 相 界 面 的 存在 ， 采 用 Latva-Kokko 和 
Rothmant 所 提出 的 算式 重新 标 色 算法 ， 该 算法 能 
够 保持 界面 的 各 问 同 性 同时 有 效 减 少 相 界 面 处 的 假 
拟 速度 号 。 为 提升 模型 的 数值 稳定 性 ， 采 用 多 松弛 
时 间 (MRT) 模 型 09 处 理 碰撞 项 。 图 1 中 计算 域 上 边 
界 采 用 半 步 长 反弹 实施 无 滑 移 边界 条 件 ， 下 边界 为 
对 称 边界 条 件 ， 左 右边 界 为 周期 性 边界 条 件 。 
根据 Melcher 和 Taylor”! LAJ Saville 等 EC 对 电 
流体 动力 学 的 研究 , 由 于 液 滴 表 面 的 动态 电流 很 小 、 
磁 效 应 的 影响 可 以 忽略 ， 因 此 电场 强度 EE 是 无 旋 的 
(Vx 五 =0)， 并 且 与 电势 的 关系 可 表述 为 
E =-V0 (4) 
流体 系统 采用 漏电 介质 模型 描述 ， 即 认为 液 滴 
具有 一 定 的 导电 能 力 ; 由 于 电荷 在 界面 上 所 聚集 的 
速度 远 高 于 流体 运动 的 时 间 尺 度 ， 所 以 可 对 电荷 控 
制 方程 做 准 静 态 人 简化， 即 
V:-(cE)=0 (5) 
将 方程 (4) 代 入 到 方程 (5), 可 获得 流体 的 电荷 守恒 方 
程 为 


V-(oV®)=0 (6) 
对 于 不 可 压缩 流体 ， 作 用 在 相 界 面 处 的 电 应 力 
Fs TI'S: 


1 | 
Fs=V.r" =-7EEVetgE (7) 


其 中 ，z” 为 麦克 斯 书 应 力 张 量 ，9" =V.(E 巨 ) 为 
电荷 在 界面 处 的 体 密度 。 在 具体 实施 过 程 中 ， 首 先 
采用 有 限 差分 法 离散 求解 方程 (6) 获 得 电势 ， 然 后 通 
过 方程 (4) 获 得 电场 强度 , 进一步 通过 方程 (7) 计 算 作 
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用 在 相 界面 处 的 电 应 力 并 与 格子 Boltzmann 颜色 模 
型 相 耦 合 。 


2 模型 验证 


本 节 通 过 模拟 电场 毛细 数 Car 对 液 滴 变 形 的 影 
响 ， 并 与 Feng 等 中 提出 的 理论 解 相 比 较 ， 验 证 所 
发 展 模型 的 准确 性 。 不 失 一 般 性 ， 考 虑 三 组 典型 的 
(R,O) 值 ， 分 别 为 (1.75，3.5)，f,= 一 4.69<0; (3.25, 
3.5), f=4.31>0; (4.75，3.5)，_ff=17.81>0。 对 于 
每 组 (R,O) 值 ， 研 究 液 滴 变形 参数 D 随 电场 毛细 数 
Cag(0< Ca, <1) 的 变化 规律 。 


0.2 上 — R=1.75 Theory 
—— R=3.25 Theory 
—— R=4.75 Theory 


| m R=1.75 Simulation 
@ R=3.25 Simulation 
-0.2; 全 R=4.75 Simulation 


0 0.2 0.4 0.6 0.8 1 


图 2 不 同 (R,O) 情 况 下 ， 液 滴 变形 参数 DD 随 电 场 毛 细 数 Cay 
的 变化 。 实 线 为 Feng 等 [5 的 理论 解 ; 符号 为 本 文 的 模拟 结 
果 


Fig.2 Deformation parameter D as a function of Caz, at different 


(R,Q) values. The solid lines are the theoretical values 
calculated from Feng et al. !'*!, while the discrete symbols are 


the simulation results 


从 图 2 可 知 ， 对 于 任何 一 组 (R,0) 值 ， 液 滴 变 形 
(扁平 或 遍 长 ) 程 度 随 电场 毛细 数 Cas 的 增加 而 增加 。 
电场 毛细 数 Cas 小 于 0.3 时 ， 模 拟 结果 与 理论 解 吻 
合 良好 ;， 然 而，Cas 大 于 0.3 时 ， 对 于 R=1.75 或 
R=4.75， 由 于 液 滴 变 形 较 大 ， 不 再 满足 理论 推导 中 
所 做 的 小 变形 假设 ， 模 拟 结 果 与 理论 解 发 生 偏差 ， 
并 且 偏 差 随 Cas 的 增 大 而 增 大 。Hu 等 人 PR" 对 于 三 
维 液 滴 变形 的 数值 模拟 也 得 出 与 之 相似 的 结论 。 图 
3 给 出 了 四 种 参数 条 件 下 液 滴 稳 态 变 形 情况 、 流 线 
和 速度 矢量 图 。 从 中 可 以 看 出 ， 扁 平和 扁 长 型 液 滴 
内 的 环流 方向 刚好 相反 。 


以 上 研究 结果 表明 ， 本 文 数值 模型 在 恒 稳 电场 
作用 下 液 滴 小 变形 的 数值 模拟 中 ， 体 现 了 良好 的 数 
值 精度 和 模型 稳定 性 。 


图 3 不 同 参数 条 件 下 稳 态 的 液 滴 形 状 、 流 线 ( 左 ) 及 速度 矢量 
分 布 ( 右 ): (a) R=1.75，O=3.5，Car=0.5; (b) R=1.75，O=3.5， 
Cag=1.0; (c) R=4.75, Q=3.5, Cag=0.5; (d) R=4.75, Q=3.5, 
Cam=1.0。 其 中 ， 蓝 色 实 线 为 液 注 相 界面 

Fig.3 The steady droplet shape, streamlines (left half) and the 
velocity vectors (right half) for (a) R=1.75, Q=3.5, Ca,=0.5; 
(b) R=1.75, Q=3.5, Cag=1.0; (c) R=4.75, Q=3.5, Cag=0.5; 

(d) R=4.75, Q=3.5, Ca,=1.0. The blue lines represent the 


droplet interfaces 
3 结果 与 讨论 


3. 1 液 滴 扁平 型 变形 及 破裂 

指示 函数 <0 时 ， 液 滴 会 发 生 扁平 型 变形 ， 随 
着 电场 毛细 数 Car 的 进一步 增 大 ， 液 滴 将 可 能 发 生 
破裂 。 图 4 给 出 了 R=0.1, O=2.0 时 的 液 滴 变形 参数 
D 随 电场 毛细 数 Cag 的 变化 情况 。 

从 图 4 可 以 看 出 ， 液 滴 变 形 参 数 DD 值 为 负数 ， 
表明 液 滴 治 着 与 电场 强度 垂直 的 方向 变形 ， 并 且 随 
着 电场 毛细 数 Cas 的 增加 而 增 大 。 根 据 毛细 数 的 大 
小 ， 我 们 观察 到 四 种 不 同 的 液 滴 形态 ， 椭 圆 形 、 类 
椭圆 形 、 哑 铃 形 和 液 滴 破裂 。 当 电场 毛细 数 
Car<0.04 时 , 稳 态 液 滴 为 椭圆 形 , 模拟 结果 与 Feng 
等 0 的 理论 解 吻合 良好 :， 当 电 场 毛 细 数 
0.04<Car<0.14 时 ， 稳 态 液 滴 为 类 椭圆 形 , 模拟 结果 
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开始 与 理论 解 出 现 差异 ( 因 变 形 较 大 , 小 变形 假设 不 
再 适用 ); 当 0.14< Cag<0.17 时 , 液 滴 长 轴 继 续 增 大 ， 
短 轴 进 一 步 减 小 ， 稳 态 液 滴 为 治 着 水 平方 向 的 哑铃 
形 。 值 得 强调 ， 当 电场 毛细 数 Cas=0.14 左右 时 ， 液 
滴 变形 参数 DD 发 生 了 突变 , 液 滴 表 面 的 电 应 力 不 再 
能 被 液 滴 内 部 环流 所 产生 的 粘性 应 力 所 平 衡 ， 液 滴 
被 拉 伸 为 两 头 大 中 间 扁 的 哑铃 型 ， 直 到 电场 力 与 流 
场 的 切 应 力 达 到 新 的 平衡 。 图 5 给 出 了 Core=0.139 
及 0.168 条 件 下 稳 态 液 滴 形 状 、 流 线 ( 左 ) 和 速度 矢量 
分 布 图 ( 右 )， 此 时 液 滴 长 宽 比 大 于 7。 从 中 可 以 看 出 
液 滴 在 水 平方 向 上 的 长 度 随 电场 强度 的 增加 而 增 
长 ， 且 液 滴 的 变形 参数 DD 趋 于 -1。 

当 电 场 毛细 数 增加 到 Cap=0.17 左右 时 , 液 滴 在 
恒 稳 电场 的 作用 下 发 生 破裂 。 如 图 6 所 示 ， 在 电场 
力 的 作用 下 ， 液 滴 中 部 发 生 收缩 ， 沿 水 平方 向 整体 
拉 伸 为 哑铃 型 。 在 液 滴 破 裂 过程 中 ， 液 滴 未 在 竖 直 
中 轴线 处 发 生 断 裂 ， 而 是 在 左右 哑铃 型 的 颈 部 发 生 
断裂 。 断 裂 后 形成 两 个 对 称 的 大 液 滴 和 中 间 一 个 小 
液 滴 ， 中 间 的 小 液 滴 会 在 界面 张力 的 主导 下 重新 收 
缩 为 椭圆 形 ， 而 左右 的 两 个 大 液 滴 会 在 电场 力 的 作 
用 下 拉 伸 为 不 对 称 的 哑铃 型 。 


Q 
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图 4 R=0.1, O=2.0 时 ， 液 滴 变 形 参 数 D 随 电场 毛细 数 Caz 
的 变化 。 黑 色 实 线 为 Feng 等 05 的 理论 解 ， 红 色 离 散 点 为 本 
文 的 模拟 结 


Fig.4 Deformation parameter D as a function of the 


electrical capillary number Ca, for R=0.1 and Q=2.0. The black 
line is the theoretical results from Feng et al. U5] while the red 


circles are the simulation results 
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(a)Caz =0.139 


(b)Cay =0.168 


图 5 R=0.1, O=2.0 时 ， 不 同 电场 毛细 数 下 稳 态 液 滴 形 
状 、 流 线 ( 左 ) 和 速度 矢量 分 布 ( 右 ); (a) Cag=0.139, (b) 
Cas=0.168。 其 中 蓝 色 实 线 为 液 滴 相 界面 
Fig.5 The steady droplet shape, streamlines (left half) and 
the velocity vectors (right half) at R=0.1 and Q=2.0 for (a) 
Ca,=0.139 and (b) Caz=0.168. The blue lines represent the 


droplet interfaces 


Q 


C > 
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图 6 Car=0.176 时 液 滴 破 裂 过 程 图 
Fig.6 The snapshots of the droplet breakup when 
Ca;=0.176 
3. 2 液 滴 扁 长 型 变形 与 破裂 

HERRA 时 ， 液 滴 会 发 生 遍 长 型 变形 。 
本 节选 取 R=10, Q=0.2 研究 不 同 电场 毛细 数 下 液 滴 
局 长 型 变形 与 破裂 行为 。 

在 模拟 中 我 们 发 现 ， 液 滴 在 电场 毛细 数 Cas 从 
零 增 大 的 过 程 中 会 呈现 四 种 不 同 的 形态 : 椭圆 形 、 
类 椭圆 形 、 液 滴 振 荡 和 液 滴 破裂 。 如 图 7 所 示 ， 电 
场 毛细 数 Cas 较 小 时 (Cas<0.2)， 液 滴 变 形 也 较 小 ， 
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沿 电场 强 度 方 向 被 拉 伸 为 椭圆 形 ， 与 理论 预测 吻合 
恨 好 。 随 着 电场 毛细 数 Car 的 增 大 , 液 滴 变 形 增 加 ， 
逐渐 偏离 彬 圆 形 (类 椭圆 形 液 滴 )。 如 图 8 所 示 ， 当 
0.6< Cag<0.88 时 ， 液 滴 变 形 参 数 D 经 历 类 阻尼 振 
荡 , 即 液 滴 在 经 过 一 段 时 间 后 最 终 趋 向 于 平衡 状态 ， 
而 且 Car 越 大 ， 液 滴 变 形 参数 万 的 峰值 越 高 ， 达 到 
稳 态 所 需 的 时 间 越 长 。 随 着 Cas 继续 增加 
(Cap>0.88)， 扁 长 型 液 滴 的 变化 与 扁平 型 液 滴 的 变 
化 出 现 显著 不 同 : 局 长 型 液 滴 在 恒 稳 电场 的 作用 下 
发 生 周 期 性 振荡 , 且 振荡 频率 与 振幅 不 随时 间 改 变 。 
该 物理 现象 在 以 往 研究 中 未 见报 道 。 当 周期 性 振荡 
出 现时 ， 随 Car 的 增 大 ， 液 滴 变 形 参 数 D 的 振荡 平 
均值 增 大 ， 振 荡 频 率 稍 许 提 高 ， 振 幅 减 小 。 


0.57 


—— Theory 
—O- Simulation 
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图 7 R=10, O=0.2 时 ， 液 滴 变 形 参 数 D 随 电场 毛细 数 
Car 的 变化 。 黑 色 实 线 为 Taylor 理论 解 ， 红 色 离 散 点 为 本 文 
模拟 结果 
Fig.7 Deformation parameter D as a function of the electrical 


capillary number Ca, for R=10 and Q=0.2. The black line is the 


theoretical results from Feng et al. "5l, while the red circles are 


the simulation results 


图 9 展示 了 电场 毛细 数 Cas=1.00 时 的 液 滴 振 荡 
变形 图 ， 其 中 图 9(a) 给 出 了 半 个 振荡 周期 内 5 个 等 
间隔 时 刻 ( 见 图 9(b)) 液 滴 的 界面 形状 。 由 图 9(a) 可 
知 ， 液 滴 呈 现 哑铃 形状 的 时 间 为 (t1-t4)， 远 长 于 呈 
现 椭圆 形状 (t5) 的 时 间 。 

和 置 于 恒 稳 电场 中 的 静止 液 滴 ， 会 在 液 滴 相 界 面 
处 产生 感应 电荷 ， 电 场 提供 给 液 滴 的 能 量 转 化 为 液 
滴 相 界 面 处 流体 的 动能 ， 使 得 液 滴 被 逐渐 拉 长 。 随 
着 液 滴 在 电场 方向 上 的 拉 伸 ， 液 滴 表 面 由 恒 稳 电场 
存在 而 产生 的 感应 电荷 量 逐 渐 减 少 ， 电 场 力也 随 之 
减弱 。 与 此 同时 ， 随 着 液 滴 变形 量 的 增加 ， 作 用 在 
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相 界 面 上 的 界面 张力 不 断 增 大 。 在 液 滴 被 拉 长 的 过 
程 中 ， 存 在 某 一 时 刻 ， 界 面 张力 对 液 滴 的 作用 超过 
电场 力 的 作用 ， 而 且 液 滴 的 动能 也 因 液 滴 内 外 环流 
的 存在 而 逐渐 耗 散 ， 动 能 耗 散 为 零 时 ， 液 滴 变形 参 
数 达 最 大 值 ( 见 图 9 中 t 时刻)。 随 后 ， 液 滴 在 界面 
张力 的 作用 下 收缩 ， 随 着 液 滴 的 收缩 ， 界 面 张力 逐 
渐 减 小 ， 液 滴 表 面 的 感应 电荷 增加 ， 电 场 力 增 大 ， 
再 次 达到 液 滴 动能 为 零 时 , 液 滴 变 形 参 数 D 达 最 小 
值 ， 液 滴 又 将 在 电场 力 的 作用 下 被 拉 长 ， 重 复 上 述 
过 程 。 因 此 ， 在 某 一 电场 毛细 数 Cas 的 范围 内 ， 液 
滴 在 电场 力 与 界面 张力 的 共同 作用 下 发 生 周期 振 
荡 。 
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图 8 液 滴 变 形 参 数 D 随时 间 步 t 的 变化 曲线 


Fig.8 The droplet deformation D versus time step t 


t x10 
(b) 
图 9 Cag=1.00 时 ，(a) 几 个 典型 时 刻 的 液 滴 形 状 ; (0) 
期 性 振荡 液 滴 的 变形 参数 随时 间 的 演化 。 其 中 (a) 图 中 各 时 刻 
在 D-Ca 曲线 上 的 位 置 已 标注 于 (b) 图 ; ty 对 应 初始 液 滴 未 发 
生 形变 的 时 刻 
Fig.9 (a) The droplet shapes at several typical times and 


(b) the time evolution of the deformation parameter D during 
the droplet oscillation when Ca;=1.00. Note that the times in (a) 
are labelled on the D-Cag curve in (b). The black lines represent 


the initial shape of the droplet 
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电场 毛细 数 Cag 增加 到 足够 大 时 (Cass>2.24)， 
液 滴 将 在 恒 稳 电场 的 作用 下 发 生 破 裂 。 对 于 扁 长 型 
变形 的 液 滴 ， 当 沿 电场 强度 方向 被 拉 长 时 ， 其 水 平 
方向 受 力 平 衡 ， 并 在 水 平 对 称 轴 上 的 液 滴 界 面 宽 度 
趋 于 零 时 发 生 破 裂 , 如 图 10 所 示 。 液 滴 破 裂 为 三 个 
子 液 滴 ， 在 中 轴线 处 有 少量 的 液 滴 残 余 ， 上 下 等 量 
的 两 部 分 被 拉 伸 为 长 条 形 并 沿 着 电场 方向 移动 。 值 
得 指出 ， 扁 长 型 液 滴 破 裂 时 的 变形 参数 DAN 
0.99( 即 液 滴 的 长 宽 比 约 为 100:1), 该 参数 ( 子 液 滴 的 
大 小 ) 几 乎 不 随 Cag 的 增加 而 变化 , 而 液 滴 破 裂 所 需 
时 间 随 Cag 的 增 大 而 减 小 。 


图 10 Cas=2.26 时 ， 液 滴 破 裂 过 程 图 
Fig.10 The snapshots of the droplet breakup when 
Cag=2.26 


基于 上 述 讨论 ， 针 对 扁 长 型 变形 的 液 滴 ， 可 以 
得 出 以 下 结论 : 在 恒定 介 电 常数 比 O 和 电导 率 比 R 
的 情况 下 ， 液 滴 存 在 两 个 电场 临界 毛细 数 ， 分 别 为 
区 分 液 滴 发 生 稳 态 变形 或 发 生 周 期 性 振荡 的 临界 毛 
细 数 Carr, 和 液 滴 破 裂 的 最 小 电场 毛细 数 Cagpp。 对 
于 当前 算 例 R=10, O=0.2), 振荡 临界 毛细 数 Cagy E 
区 间 (0.87,0.89) 中 ， 破 裂 临界 毛细 数 Cars 在 区 间 
(2.24,2.26) 中 。 


4 结论 
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本 文 基 于 格子 Boltzmann 颜色 模型 计算 两 相 流 
动 、 有 限 差 分 方法 计算 电场 ， 并 通过 漏电 介质 模型 
合 ， 发 展 出 混合 数值 方法 对 恒 稳 电场 作用 下 液 滴 
变形 与 破裂 行为 进行 研究 。 针 对 不 同 拉 伸 方向 的 液 
滴 变 形 数值 研究 发 现 ， 液 滴 变 形 (扁平 或 扇 长 ) 均 随 
电场 毛细 数 的 增加 而 增加 ; 电场 毛细 数 较 小 时 
(Cae<0.3)， 本 文 模拟 结果 与 Feng 等 03 的 理论 解 完 
全 吻合 ， 电 场 毛 细 数 较 大 时 ， 由 于 液 滴 变 形 不 符合 
小 变形 假设 ， 模 型 结果 与 理论 解 出 现 偏差 ， 并 且 偏 
差 随 Cap 的 增 大 而 增 大 。 

本 文 也 研究 了 电场 毛细 数 对 液 滴 变 形 及 破裂 现 
象 的 影响 规律 。 对 于 扁平 型 变形 的 液 滴 (f<0, R=0.1, 
2=2.0)， 液 滴 沿 着 垂直 于 电场 强度 的 方向 上 拉 伸 ， 
液 滴 随 电场 毛细 数 的 增加 呈现 四 种 不 同 的 形态 : ABE 
圆 形 (Cas<0.04)、 类 椭圆 形 (0.04<Cas<0.14)、 哑 铃 
型 (0.14< Cas<0.17) 及 破裂 (破裂 临界 毛细 数 Cap 约 
为 0.17)。 对 于 局 长 型 变形 的 液 滴 (f20，R=10， 
Q=0.2)， 液 滴 沿 着 电场 强度 方向 发 生 拉 伸 变形 ， 其 
变形 参数 万 随 电 场 毛细 数 的 变化 与 扁平 型 液 滴 显著 
不 同 ， 液 滴 呈 现 的 形态 有 : 椭圆 形 (Cas&<0.2)、 类 李 
圆 形 (0.2<Cas< 0.88). Jal Hye (0.88<Cags 2.24) 
破裂 (Cags>2.24)。 液 滴 变 形 阶 段 ， 扁 平 或 扁 长 型 液 
滴 变形 参数 D 均 随 电场 毛细 数 的 增加 而 增 大 ; 液 滴 
破裂 阶段 ， 破 裂 所 需 时 间 均 随 电场 毛细 数 的 增加 而 
缩短 。 
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